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Abstract 

Stable non-top ological solitons, Q-balls, are studied using analytical and nu- 
merical methods. Three different physically interesting potentials that support 
Q-ball solutions are considered: two typical polynomial potentials and a loga- 
rithmic potential inspired by supersymmetry. It is shown that Q-balls in these 
potentials exhibit different properties in the thick-wall limit where the charge of 
a Q-ball is typically considerably smaller than in the thin- wall limit. Analytical 
criteria are derived to check whether stable Q-balls exists in the thick-wall limit 
for typical potentials. Q-ball charge, energy and profiles are presented for each 
potential studied. Evaporation rates are calculated in the perfect thin-wall limit 
and for realistic Q-ball profiles. It is shown that in each case the evaporation 
rate increases with decreasing charge. 
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1 Introduction 



A scalar field theory with a spontaneously broken [/(l)-symmetry may contain stable 
non-topological solitons |l| , Q-balls . A Q-ball is a coherent state of a complex scalar 
field that carries a global U{1) charge. (One should be careful not to use the word 
soliton by itself to describe these objects since Q-balls are not solitons in the strict 
sense of the word). In the sector of fixed charge the Q-ball solution is the ground 
state so that its stability and existence are due to the conservation of the U{1) charge. 
Q-balls associated with a local U{1) charge have also been constructed 0] but these 
will not be discussed in this paper. 

In realistic theories Q-balls are generally allowed in supersymmetric generalizations 
of the standard model with fiat directions in their scalar potentials. Q-balls can be 
formed from scalars carrying a conserved U{1) charge. Flat directions in the Minimally 
Supersymmetric Standard Model (MSSM) have been shown to exist and they have 
been classified in In particular Q-balls have been shown to be present in the MSSM 
[§, 0] where leptonic and baryonic balls may exist. For a Q-ball solution to be possible 
in the scalar sector of the theory with potential U{(f)), U{(f))/(f)'^ should have a global 
minimum at a non zero value of 0. Such renormalizable F- and D-fiat directions 
where this is possible in the MSSM are for example the if^L-direction (where the 
expectation values iJ^ and ul are non- zero and equal) and the w'^d^rf'^-direction (along 
which Uc, dc and d'^ VEVs are non-zero) 0. For a large a D-fiat direction in the 
potential (suppressed at the Planck scale) is of the form [0 



,2 , A^l^pC-^-i) , ^AxXr 



U{ct>) = mi,\<p\^ + + + h.c. ), (1) 

Mpi ^yJ-pi 

where ms is the soft SUSY breaking mass parameter, A\ the A-term and d is the 
dimension of the non-renormalizable term in the superpotential. For Q-balls to exist 
in the theory, ms must be dependent and furthermore be such that f/(0)/|0p has a 
global minimum at a non- zero 0. 

Q-balls may be formed in the early universe by a mechanism |^, |^ that is closely 
related to the Affleck-Dine baryogenesis |^ . A scalar condensate first forms along a fiat 
direction in the potential. In the MSSM this is a D-flat direction composed of squark 
and possibly slepton fields. The spatially homogeneous condensate then develops an 
instability that leads to the formation of areas of different charge densities. This charge 
pattern can then lead to the creation of Q-balls (recent work [l^ shows that initially 
the objects formed from the AD-condensate are not in general Q-balls, but Q-axitons). 
In the MSSM Q-balls formed in such a way are baryon (and lepton) number carrying 
balls with charges typically larger than 10^'' [0. The subsequent evolution of the Q- 
balls then critically depends upon the mechanism how supersymmetry is broken in 
the theory ||lT|. If SUSY breaking occurs via a gauge mediated mechanism [|12| the 



potential for the scalar condensate is completely flat {ms behaves as for large 



enough in this scenario) and the resulting Q-balls can have huge charges |1^ . For 
a large enough Q-ball with no lepton number, all decay modes can be excluded and the 
Q-ball can be completely stable. If, however, the SUSY breaking mechanism is via a 
supergravity hidden sector [Q the scalar potential is not completely flat but radiative 
corrections allow for Q-ball solutions to exist in the potential. The mass parameter 
is in this case given by [0 

ml^ml{l + K\og{^-^)), (2) 

where K is typically —0.01 - —0.1 and Mx is a large mass scale (m|(Mx) = mo). 
In this case the Q-balls are generally unstable and can typically decay to quarks and 
nucleons [|Tl] . 



The possible existence of Q-balls has inspired new ideas and particularly their cos- 
mological significance has been discussed recently. The question of Q-ball stability 
plays an important role in their importance to cosmology. If stable Q-balls are formed 
in the early universe they can contribute to the dark matter content of the universe 
[§, 01 . These can be huge balls with charges of the order of 10^° |^ but also very 



small Q-balls can be considered as dark matter WM. Stable Q-balls can also catalyze 
explosions of neutron stars by accumulating inside the stars and decreasing the mass 
of the star by absorbing baryons ll5|. Decaying Q-balls can also be of crucial cos- 
mological significance. If Q-balls decay after the electroweak phase transition, they 
can protect baryons from the erasure of baryon number due to sphaleron transitions 
PH . Possibly they can also create the asymmetry from a. B — L conserving condensate 
[0]. Furthermore the Q-ball decay may result in the production of dark matter in the 
form of the lightest supersymmetric particle (LSPs). This process can then naturally 
explain the baryon to dark matter ratio of the universe |ll|]. Small Q-balls produced 



in a collider would also have interesting experimental properties ||14|. It has even been 
suggested that Q-ball catalyzed processes may provide us with an inexhaustible source 
of energy |13|. 

The detection of Q-balls was discussed in |]T6[ where detection was considered 



for two types of Q-balls, Supersymmetric Electrically Neutral Solitons (SENS) and 
Supersymmetric Electrically Charged Solitons (SECS). The mode of detection depends 
crucially on the possible electric charge carried by the Q-ball. Assuming that all or 
most of dark matter is in Q-balls, the number density of stable Q-balls was estimated 



to be 16 



nQ ~ 5 X 10-'Q-'/\^^) cm-^ (3) 

where Qb is the charge of the baryonic ball and m is the mass of the baryon number 
carrying scalar. Further assuming that the average velocity of Q-balls is ~ 10~^c, 



experimental limits can then be calculated. From the MACRO search a lower 
bound on SECS charge is ~ 10^^. From monopole searches |I18| the lower bound for 



SENS charge is ~ 10^^. In considering these limits one must bear in mind that here 
it has been assumed that dark matter consist mostly of Q-balls. If there are far fewer 
Q-balls in the universe, however, then this bounds can become significantly lower. 
Also if Q-balls are unstable and evaporate there may be only indirect evidence of their 
existence present in the universe today. 

In this paper some of the analytical and numerical properties of Q-balls are studied. 
We first derive some simple analytical results for Q-balls in generic scalar potentials. 
Q-balls are then studied numerically in these typical potentials. In the last part of this 
paper we consider Q-ball evaporation for realistic Q-ball profiles and present numerical 
results. 

2 Q-ball solutions 

We consider a field theory with a f/(l) symmetric scalar potential. U{(j)) has a global 
minimum U{0) = at = and the scalar field has a unit charge with respect to 
the f/(l). The Lagrangian density for the theory is given by 

/: = a^0*9>-f/(0). (4) 

The charge and energy of a given field configuration are 

Q = h I (0*5*0 - <pdt(t)*)d\ (5) 
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E= I d'x[^ 0- + l|V0|Vt/(0)]. (6) 
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We wish to find the minimum energy configuration at a fixed charge, i. e.the Q-ball 
solution. If it is energetically favourable to store charge in a Q-ball compared to free 
particles, the Q-ball will be stable (in realistic theories, however, additional couplings 
may lead to evaporation in a different sector of the theory). Hence for a stable Q-ball, 
condition E < mQ, where m is the mass of the 0-scalar, must hold. Minimizing the 
energy is done by using Lagrange multipliers i.e.we wish to minimize 

e^ = E + u[Q-^J {ct>*dt(t> - 09t0*)rf3x], (7) 

where is a Lagrange multiplier. Separating (|^) into time and space dependent parts 
we get 

t^ = j d'x]^ \dtct> - iujct>f + j d'x[]^ I V0I' + UM)] + ^Q, (8) 



where 

UM = U{<P) - \uJ^\ (9) 
To minimize (|^) we first choose 

0(x,t) = e'-V(a;), (10) 

where 0(a;) is now time independent and real. Charge and energy are now given by 

Q = ujJ (t){xf<fx (11) 



and 



E= / rf=^x[^|V0(x)|' + t/^(0(x))]+a;g. (12) 



The problem of minimizing the energy for a fixed uj in four dimensions is then equiv- 
alent to the problem of finding the bounce for tunneling in D = 3 If a bounce 



solution exists it is spherically symmetric [0. Hence the Q-ball solution also has a 
spherical symmetry, 

0(x) = 0(r). (13) 

The charge and energy expressions for a spherically symmetric Q-ball are given by 

Q = Atiu r^(f){rfdr (14) 
Jo 

and 

E = 47r []^{^f + UM{rWdr + ujQ. (15) 
The equation of motion at a fixed uj is 

+ (16) 

To obtain the Q-ball profiles we must solve (p!6D with boundary conditions 0'(O) = 
0, (t){oo) = 0. The initial value of the field, 0(0), is left as a free parameter. This is, 
in general, a numerical problem that is discussed later in this paper. 

Some properties of Q-ball solutions in general potentials can be studied analytically. 
These are discussed in the following. 

2.1 Behaviour of Q-ball charge as u ^ m 

We examine a general potential of the following form in D dimensions, 

Ui<P) = \{m' - u?)<^^ - + O(0^+^), (17) 



where A > 2 and A is positive. When r is large the 'friction' term ^cp' is small so that 
the equation of motion can be approximated by 

<P"^d^UM- (18) 

As r — > oo, tends to zero and we can neglect all but the effective mass term in the 
potential Uuj{4>) so that from (|l^) we get 

0(r) ~ e-^™'-^'". (19) 

Furthermore, since when u m the position of the zero of Uuj{(l)) must be near the 
origin [|l|, we can improve the approximation (|19D to 

0~ (m2-cu2)A^e-^'"'-"''-. (20) 

We can now calculate the charge at this limit, 

Q ~ lim / d'^xcP' ~ (m^ - u')^^ H drr^'-'e'^^^^^^ ~ (m^ - a;2)(^-f ). (21) 

ui-*mj Jo 

Clearly the value of 77 = (^^ — ^) now determines the behaviour of the Q-ball charge 
as — m; 

> g ^ 

= ^ Q ^ const. (22) 
?7 < =^ Q ^ 00. 

This result agrees with the results for A = 4 presented in 

Even though solving the Q-ball equation of motion is generally a numerical prob- 
lem, in certain cases analytical approximations can be useful. Two limiting cases are 
discussed in the following, the thin-wall and thick-wall approximations. 



2.2 Thin- wall approximation 

In a flat potential a Q-ball may have a large charge and the effect of its surface can 
be ignored compared to the bulk of the ball. Such a thin-walled Q-ball [Q may be 
approximated by an ansatz 

(23) 

so that all surface effects are ignored. 0o is defined such that U (0) / 0^ is minimum for 
= 00- The charge is given by 

4:71 

Q = = Vu<Pl (24) 




and the energy is 

E = Q 



\ 



(25) 



In a flat potential ratio E/Q then decreases for increasing (p so that it becomes more 
and more energetically favourable to store quanta in the Q-ball. This approximation 
is accurate for very large, thin-walled Q-balls. For the types of potentials studied in 
this paper, a thin-wall approximation is valid for small u and it generally deteriorates 
as uj tends to m. Generally the existence of a thin-walled Q-balls is associated with 
completely fiat directions in the scalar potential of the theory. These can lead to the 
formation of Q-balls with very large charges M, . 



2.3 Thick- wall approximation 

As UJ tends to m, the Q-ball configuration becomes more and more thick- walled so 
that at some point the thin-wall approximation is no longer valid. When u is close to 
m another approximation, the thick-wall approximation, can be used. 
Kusenko considered a thick-wall approximation in the potential 

U{<P) = ImV' - + A>^ (26) 

so that 

f/^(0) = ^{m' - u;')<j>' - \<f + Ay. (27) 
We generalize the potential to 

[/(0) = imV-A0^ + O(0^+^), (28) 

where A > 2 and A > and consider the possibility of Q-balls in D dimensions. In 
the thick-wall limit, ie.when the higher terms are negligible, (9(0^"^^) ~ 0, 

U^(ct>)^^{m'-u;')4>'-X<P^. (29) 

Introducing dimensionless space-time coordinates, = (M^ — uj'^)^Xi, and a dimen- 
sionless field, ip = {.j^Tz;ji)^~^'t'^ the expression for energy (|12|) can be written as 

2 A n r, 2A-AD+2D 

e^ = \—S^{m'-oj^) 2{A-2) +^Q^ (30) 

where is the action of the bounce in the dimensionless potential ^ip"^ — ip^. The 
exact value of the action is not essential to our argument provided that it is well 
behaved as D and A grow. Sp^ is plotted in Fig. for D = 3. We have also checked 
that SpL is well behaved for higher dimensions. 
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Figure 1: The three dimensional bounce for different values of A 

Q-balls may exist provided that has a minimum for some < < m such that 
^LUc < fnQ. Furthermore for the thick wall-approximation to be valid Uc must be close 
to m. 

The expression for energy (|30|) and its derivative 



Q-\—2uj{m^-u^fS^{C + l), 



(31) 



where C = "^"2(^-2)^ ' exhibit different behaviour as u; — > m, depending on the value of 



If C > 0, mQ and ^' 



provided that the condition 

Q < 2A^m2^+i(2C)^(l + 2C)-^-^^^(C + 1 

holds. 



Q. In this case Q-balls exist in the thick- wall limit 

(32) 



If = no local minimum exists in the interval < < m and hence no stable 
Q-balls exist in the thick- wall limit. 

If— 1<C<0, e^— > mQ but now ^ — > — 00 so that there are no stable Q-ball 
solutions in the thick- wall limit. 

If = — 1 no local minimum exists and hence there are no stable Q-ball solutions 
in the thick- wall limit. 



If C < — 1, ^ oo and ^ oo. Again, there are no stable Q-balls in the 
thick- wall limit. 

To summarize, for stable Q-balls to exist in the thick-wall limit, conditions 
0<A + 2D-AD 

Q < 2A^m2f+i(2C)^(l + 2C)-^-55^(C + 1), ^^^^ 

must hold. 

For example, when D = 3, A must be less than y . In it was shown that Q-balls 
exist at the thick- wall limit when D = 3, A = 3. When D = 3, A = 4 one can easily 
show that 



E = mg^l + /3-2, (34) 

where /3 = 4^. Since (3 > 0, E > mQ and hence there are no stable Q-ball solutions. 
This is also demonstrated numerically in the following section. 



3 Numerical results 

We have examined three different types of potentials numerically, namely 

UM = ^ml<p^ - a^<P^ + X^<|)^ (35) 

UM = - y 0^ + (36) 

Usi<P) = m^</)2(l - inog(^)) + A30^°. (37) 

This choice of potentials is motivated by physical arguments. Ui (0) is a typical simple 
potential in field theories that can contain extended objects of the Q-ball type. U2{(p) is 
a typical potential in finite temperature field theories as the the presence of the 0^ term 
indicates. U^{(f)) corresponds to a D-flat direction in the MSSM with SUSY broken via a 
conventional supergravity hidden sector. The parameters are chosen for each potential 
such that when = there is a degenerate global minimum, f/(0 = 0) = and there 
exists a value 0o > such that f/(0o)=O, which sets ai = y^2mf Ai, 02 = \J ^1^2X2 
and M = {jKm^X^^ exp(— 1 — ■^))^. This choice ensures that there is no spontaneous 
breakdown of f/(l)-symmetry at = 0. Furthermore the parameters must chosen 
appropriately to ensure the existence of Q-balls as u increases ^e.the coefficient of the 
highest power should be small enough compared to m. We also need to ensure the 
correct behaviour at large i.e.U{(f)) — 00 as — > 00. 

The parameter values we have chosen to plot the figures are: 

mi = 1 7712 = 1 ^3 = 100 GeV 

Ai = 0.001 A2 = 0.01 A3 = ^ (38) 
K = 0.1. 



Note that only in potential U3{(j)) we need to introduce dimensional units due to the 
presence of the Planck mass parameter. The choice of parameters is not special for 
Ui{(f)) and U2{4>), the value of m simply sets the scale of the units, Ai and A2 set the 
profile of the potential. The parameter values for U^{(j)) are typical values for a realistic 
potential. Solving for the place of the local minimum 0i and maximum 02? we get 



^(2) _ _ /2 



^(1) _ 1^(1) a(2) _ 1^(2) ^^^> 

1 



where superscripts (1, 2) refer to different potentials ?7i,2(0)- From these the height of 
the potential barrier can be calculated f/(02^'^'*) = f^. The effect of changing A 



27V2A' 16A- 

in f/i and U2 can be seen by considering the associated Lagrangian with a potential 
of the form 




U{<P)=h<P\<P^-\—)\ (40) 



where is a constant. Clearly this potential is of the same form as Ui and U2 with 
the assumption that at = there is a degenerate minimum. Under the scaling 
s(f), b s~'^^b the Lagrangian transforms as £ — s^£. Charge and energy scale 
as Q ^ IP'Q E — s> IP'E while the radius of the Q-ball is invariant as can be seen from 
the equation of motion (|T^) . From this we can see the effect of changing Aj A^ 
(i = 1, 2) in Ui and U2 while keeping m constant; charge and energy scale as 

E i^y/^E, (41) 

where A = 2,1 for f/i,2(0), respectively. The plotted figures therefore represent typical 
Q-ball solutions, changing A simply scales the field 0. 

We will now consider some general properties of Q-ball solutions. As u increases 
from zero, a new global minimum 0o with U (0o) < develops. This does not, however, 
automatically imply the existence of Q-balls as can be seen from the equation of 
motion ([T6|). If we consider the problem of solving the equation of motion in terms of 
a mechanical analogy, where a Q-ball solution is analogous to the motion of a particle 
in the potential — f/(0) with r playing the role of time, we can easily see that the 
presence of the 'friction' term requires that u must be large enough to ensure that 
the initial 'potential energy' of the particle is large enough for the particle to reach 
= 0. Put more mathematically, integrating (|T6|) with respect to r and using the fact 
that U^{(f){oo)) = one gets 

roo 1 

[/^(0(O)) + 2/ rfr-(0'(r))2 = 0, (42) 
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Figure 2: as a function of r, U = Ui{(f)) (mi = 1) 



from which it is clear that U^{(f){0)) must be large enough to overcome the contribution 
from 'friction'. As u tends to m the potential barrier diminishes until at = m it 
disappears. As the barrier disappears Q-ball solutions no longer exist since the field 
now cannot stop at = but will continue rolling and acquire a negative value. 

Numerically the solution can be found by taking into account the boundary con- 
ditions 0'(O) = and 0(oo) = 0. The initial value of the field, 0(0), is then varied 
until a satisfactory solution is found. Depending on the value of u and the form of the 
potential, different types solutions can be found. The solutions vary from a thin-walled 
Q-ball to a thick-walled type as u tends to m in each case. This is as expected since 
as u! starts to increase a global minimum develops and 0(0) lies near the minimum 
leading to a thin-walled solution. As u increases further 0(0) moves further away from 
the global minimum leading to thick-walled Q-balls. 

We have plotted the solutions for different values of u) for each potential for a 
typical choice of parameters in Figures |^, ^ and From the Figures |^ and § the 
transformation from a thin-walled ball to a thick-walled one can clearly be seen. 

The charge and energy of a Q-ball are given by ([I^ , p!5|) . We have plotted the charge 
and energy versus u, E as a function of Q and E/Q versus u for the different types of 
potentials in Figures |5|-p!3[ In the figures where E is plotted versus Q, Figures 0, p!0| 



and |T3|, a straight line E = niQ indicating stability is drawn, Q-balls are stable (with 
respect to decays into scalars) under this line. In all the cases there is a different 
behaviour as uj tends to m. 



20 



15 



10 








co=0.3 
co=0.4 
(0=0.5 
00=0.6 
co=0.7 
(0=0.8 
(0=0.9 







10 



15 



200 



Figure 3: as a function oir, U — C/2(0) ("^2 = 1) 
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Figure 4: as a function oir, U = U^{(1)) 
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Figure 5: Energy and charge versus u, U = Ui{(f>) (mi = 1) 

For Ui{(f)) there exists a region of small u where stable thin-walled Q-balls exist. 
As uj increases charge decreases and Q-balls become thick-walled. As charge continues 
to decrease Q-balls first become unstable. Charge and energy then reach a minimum 
and start to grow again, while still in the region of instability. As u tends to one, 
E tends asymptotically to mQ. This behaviour is also clear from the E/Q versus u 
figure. Fig. |^, where a horizontal line ai E/Q = m represents the stability line. 

Even though the Q-ball profiles in the second case, U = U2{4>), look quite similar 
to those of the previous case, the charge and energy have a different behaviour as u 
increases. The profile still mutates from a thin-walled ball to a thick-walled one like 



before but now the energy versus charge plot. Fig. |10|, shows a new type of behaviour. 
Clearly E/Q < m for all u indicating that Q-balls are stable in this potential. The 
difference in Q-ball stability at the thick- wall limit in the two potentials f/i(0) and 
U2{(f)) is in accordance with the previous analytical reasoning (pSl) . 

The third case is, again, different from the two previous ones. The stability curve. 



Fig. [T^, is now such that as u increases, Q-balls become unstable like in the first case 
but here the charge and energy simply decrease with increasing u for the whole range 
of parameter values. 
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Figure 6: Energy-charge ratio versus u, U — C/i(0) (mi = 1) 
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Figure 7: Energy versus charge, U = Ui{(j)) (mi = 1) 
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Figure 8: Energy and charge versus uj,U — C/2(0) (m2 = 1) 
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Figure 9: Energy-charge ratio versus cu, U = t/2(0) (^2 = 1) 
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Figure 10: Energy versus charge, C/ = C/2(0) (m2 = 1) 



> 2 
CD 



o 







> 
CD 



.9^ Q 



CD 



-2 



OJ 

o 



-4 



30 50 70 90 

(o/GeV 



Figure 11: Energy and charge versus uj, U — U^^ct)) 
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Figure 12: Energy-charge ratio versus a;, C/ — Us{(l)) 
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Figure 13: Energy versus charge, U = U^^cj)) 



The numerical studies of these three types of potentials show different behaviour of 
charge and energy with changing uj for each type. As analytical results in the previous 
section showed the behaviour in the thick-wall limit depends critically on the power 
of the next to leading term. The logarithmic potential also shows a different kind 
of behaviour in the thick- wall limit. The question of stability of Q-balls in the thick 
wall-limit is hence strongly dependent on the exact form of the potential. 



4 Q-ball evaporation 

In realistic theories where the 0-field couples to other fields the simple criterion for 
absolute stability E < niQ no longer holds (it is still obviously true for (p decays). 
Here we discuss a decay mode to massless fermions. The evaporation of Q-balls to 
neutrinos was discussed in detail in It was shown that for a large Q-ball with a 



step-function boundary there exists an upper bound of the evaporation rate per unit 
area. 



We follow closely the approach described in ||2T[. This is a leading order semi 
classical approximation where massless fermions are considered in the presence of a 
classical field, the Q-ball background. Since the Dirac sea is filled inside the Q-ball, 
evaporation can only progress through the Q-ball surface and is hence proportional to 
the area of the Q-ball. Evaporation progresses through pair-production at the Q-ball 
surface so that each pair produced decreases the charge of the Q-ball by two (e.g. a 
lepton number carrying ball could evaporate into neutrinos so that AL = 2 for each 
pair of neutrinos produced). 

The Lagrangian density is given by 

C = d^(f)*dy - f/(|0|) + ip^ido + ia ■ - ig(j)i/j^a2ilj* + ig(j)*ilj'^a2ilj, (43) 

where ip is a two component Weyl spinor corresponding to the fermion field, is the 
Q-ball scalar field and a are the Pauli matrices. The equations of motion are 

{ido + ia ■ V)ip - g(j)x = 

{tdo-icT-V)x-g(l>*^P = 0, (44) 

where x = 'icr2'^* is a Weyl spinor. As in |^ we choose normal modes such that 
ijj ~ ^hif+^'M and x ~ e'"*-"^"'"'^'-'**]. With this choice the field equations become 

(| + a;' + za- V)V^-(70X = 
(-|+cu'-za- V)x-#^ = 0. (45) 



The eigensolutions to these equations can be found using spherical spinors. We make 
the ansatz 

= gi{r)Qjirn + ifi{r)QjUm 
X = g2{r)^jirn + if2ir)nji>m, (46) 
where fljim, ^ji'm are spherical spinors 



^j,j~l/2,m — 

^j,j+i/2,m = fjT^^ ' (47) 

^ V 2j+2 ^i,m+l/2 / 

Yim are spherical harmonics, j and m are the angular momentum quantum numbers 
and /' = 2j — I and j = / ± |. 

Substituting the ansatz into the field equations ( ^3]) the radial equations are easily 
found (we set g = 1): 

('^ + u;')g,{r) + f[{r)-Mr)^ = <j>{r)g2{r) 
2 r 

{^ + uj')h{T)-g[{T)-g,{T)'^^ = 0(r)/2(r) (48) 
z r 

(-^+co')g,{r) + f^{r) + Mr)^ = <P{r)g,{r) 

(-^+cu')/2(r) + ^?^(r) + <72(r)i^ = 0(r)/i(r). 
2 r 

We must now solve these equations with appropriate boundary conditions. Near the 

origin r = we can approximate the Q-ball potential with a constant, 0(r) = 0o- In a 

suitable basis the radial equations decouple and an analytical solution can be found: 

gl{r) = Ai{7] + a)r]iji_{r]ir) + A2{r] - a)r]2ji_i'n2r) 

/{'(r) = Ai{uj + a)r]iji^{7]ir) + A2{uj - a)ji^{r]2r) (49) 

5,0(r) = Aig(t)oriiii_{riir) + A2g(t)ori2ii_{ri2r) 

f^(r) = Aig(l)or]iji^{riir) + A2g(f)omji+iV2r), 



where l± = j ± ^, a = yuj'^ ~ g'^'Po a^id 771^2 = y ± a. Here the divergent part 
(at the origin) of the general solutions has been omitted to keep the wave functions 
normalizable at r = 0. Far away from the Q-ball surface = and the equations of 
motion decouple into two pairs. An analytical solution can again be found: 

gr{r) = A,hhf}{hr) + BiKhf}{hr) 

frir) = A,hhf}{hr) + Bikihf^{hr), (50) 



where i = 1,2, h^'^^ are the spherical Hankel functions and ki^2 = f i From these 
we can identify the in- and out-moving waves since ~ e''' and /i^^^ ~ e"^''. Solving 
for the coefficients Ai and Bi we get 



1 gr{r)h[fihr) - fr{r)h^u{kir) 
h hf-}{hr)hf^{hr) - hf}{hr)hf}{hr) 



^ ^ 1 9r{r)h^S^hr)-fr{r)h^'}{hr) 

The evaporation rate can be calculated assuming that there is no incoming x 
wave since then all the reflected x Aux must have been transmutated. From the 



transmutation coefficient T the evaporation rate can be calculated 



'i-i:rf^(^i^m(K,)?. (52) 



Here it is worth noting another result obtained in |^ 



\T{k+,j)\ = \T{k.,j)\<l. (53) 

Assuming that the potential f/(0) scales as U{s(f)) = s'^U{(j)), we note that under 
the transformation — s0, ip sip, g s^^g the Lagrangian (p3D transforms as 
£ s^C so that the evaporation rate is invariant. Assuming as before a potential of 
the form (^) we note that the evaporation rate is invariant if b scales as 6 ^ s~'^'^b. 



In the step- function boundary was considered in the large R limit. The calcu- 
lation for any R is easily done by matching the in and out solutions (|^,^). Assuming 
that there is no incoming x flux, all the reflected x waves must have then been trans- 
mutated. Identifying the coefficients Ai, Bi with the incoming, reflected and transmu- 
tated waves, the transmutation rate can be calculated in a straightforward way. The 
resulting expression is quite lengthy and is not presented here. 

The maximum evaporation rate per unit area in the i? — oo limit was calculated 
in [21] by trading the angular momentum variable j for a linear momentum variable 



and using relation (|5BD. The maximum evaporation rate in this limit is given by 

, dQ _ uj^ 
^dAdr''"'"'' ~ 1927r2' 



(54) 



We have plotted the evaporation rates per unit area as a fraction of {j^)max for a 



number of i?'s in Fig. 14 



Obviously the evaporation rate is strongly dependent on R and it approaches the 



limiting proflle given in [B^ . From the Fig. IT^ it can be seen that the evaporation rate 



per unit area flrst increases with R for all g4>o/uj until it reaches a maximum. Then 




Figure 14: Evaporation rates for a step-function boundary 



the maximum value slowly starts to decrease and shift towards smaller g(j)o/uj while 
a "bump" develops around log{g(j)o/ij) ~ —0.5. The different types of behaviour at 
small and large g(f)o/uj is due to the presence of the exp{i^Juj''^ — g'^(f)QR)-teTm in the 
matching process of the wave functions in- and outside the Q-ball. When u' > gcf)^ the 
wave-function consists of oscillatory functions compared to uj' < g(f)Q when exponen- 
tially growing and decaying parts are present. As R increases this effect becomes more 
pronounced and the evaporation rate tends to a constant more quickly as g4>o/uj in- 
crease. The slightly rough features at small g4>o/^ are also due to the strong oscillatory 
R dependence. 

In a realistic step-function is not always a good approximation. It ap- 

proximates a thin-walled Q-ball well but for a thick-walled Q-ball its accuracy may 
be questioned. The profile a thick-walled Q-ball is significantly different from a step- 
function and furthermore its radius is not well defined. To model a thick-walled Q-ball 
with a step-function Q-ball one could for example choose the charges to be equal. This 
would still leave the radius and the value of the field inside the ball arbitrary. One 
should then define the Q-ball boundary to obtain a radius and hence obtain 0o- The 
accuracy of such an approximation can be questioned especially in the case of very 
thick-walled Q-balls. 

We have examined the evaporation rates numerically for realistic Q-balls profiles 
(in addition to checking that the purely numerical calculation recovers the results for a 
step-function boundary). This is done by solving the equations (|48|) numerically for a 
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Figure 15: Evaporation rate ^ versus logQ, U = Ui{(j)) (mi = 1) 

given profile. The profiles have been calculated from the different types of potentials, 
2 = 1,2,3. 

The numerical solution is found using the fact that asymptotically the Q-ball profile 
at the origin is fiat and can be approximated with a constant. The analytical solution 
is then given by (^Q]). Far away from the Q-ball surface the potential is zero and the 
solution there is given by (|50|) . The problem then reduces to finding the right initial 
conditions to obtain a solution from which we can deduce the transmutation coefficient 
i.e.we must choose Ai and A2 in (^) such that far away from the Q-ball there is no 
incoming x-wave. The coefficients for the incoming and reflected ip-wave and for the 
transmutated x-wave can then be found by using (|5T|). 

We have calculated the evaporation rate (p2D for a number of uj:s in the three 



different potentials. These are presented in Figures |T^, |1^ and as a function of 
charge. Note that here the evaporation rate is not given as a fraction of the maximum 
rate and is the total evaporation rate (instead of the evaporation rate per unit area). 
We have only plotted the evaporation rates for a parameter range where Q-balls are 
stable with respect to decay into 0-scalars. 

From the figures we see that at large Q the evaporation rate decreases with increas- 
ing Q. This holds for the whole of the studied parameter range for all the potentials so 
that a Q-ball will evaporate faster as it decreases in size. The calculation is quite time 
consuming and requires high numerical accuracy in choosing Ai and A2 correctly. This 
fact is evident in the slight non-smoothness present in the figures due to numerical 
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Figure 16: Evaporation rate ^ versus logQ, U — U2{(f)) {mi — 1) 
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Figure 17: Evaporation rate ^ versus logQ, U — C/3(0) 
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inaccuracies during the computation. 

We have fitted a curve ^ = a{\ogQ)^ + 7 through the points in each case. These 
fits are presented in the corresponding figures. 

It is cosmologically interesting to estimate the evaporation times, especially for 
the realistic potential U^{(j)). Let us consider a typical Q-ball in the logarithmic po- 
tential that is stable with respect to scalar decays. For a Q-ball to survive past 
the electroweak (EW) phase transition we can estimate that its lifetime must be 
~ 10~^^ s ~ 10^^ (GeV)~^. From the Fig. |l^ we can estimate that the evapora- 
tion rate is ~ 10 GeV so that its charge must be at least 10^^ for it to evaporate for 
10~^^ s. However, from Fig. [1^ we see that scalar decays begin at about Q ~ 10^^. 
The lifetime of a typical Q-ball with charge of the order of ~ (M/(GeV))^ is then 
sufficiently long for it to survive past the EW phase transition. On the other hand a 
Q-ball with charge ~ 10^^ will evaporate in ~ 10^^ (GeV)~^ ~ 10~^ s so that Q-balls 
are no longer present in the universe. 

We can now estimate the emitted power in Q-ball evaporation. Writing 

dt dQ dt ^ ^ 

so that we see that |^ must be considered as E decreases. From the Figures |^, |10| and 
TB| we can see that in the stable regime, |^ < m in all the three cases. The emitted 



power then increases until scalar decays begin. Especially in the case of potential 
U2{(p) the emitted power can be quite large. However, it must be kept in mind that 
this calculation is based on a semi classical analysis and very small Q-balls should be 
analyzed quantum mechanically. 

In analyzing the decay of Q-balls in the early universe one should also consider dis- 
sociation and dissolution |Tl|. These are due to the presence of a thermal background 
in the early universe. To evade the dissociation and dissolution of typical Q-balls 
produced from an AD-condensate, an upper bound for the reheating temperature of 
~ 10'^"^ GeV has been derived UTT . 



5 Conclusions 

Q-ball properties have been studied using both analytical and numerical methods. 
Analytically we have examined the stability and charge of Q-balls in the limit where 
ui ^ m. In this thick wall limit we have derived conditions for stable Q-balls to exist 
in a potential that can be expressed as powers of in the limit of small 0. It was 
shown that stability depends on the dimension of the theory as well as on the power 
of the next to leading term in the small (p limit. These results confirm previous results 
and present a method of analyzing the stability of Q-balls when uj m. 



Three different types of potentials have been studied numerically. These calcula- 
tions show how stability in the thick-wall limit is strongly dependent on the exact form 
of the scalar potential. Numerical examinations show different types of behaviour in 
the charge versus energy plot for all of the three potentials. Furthermore, numeri- 
cal work agrees with the analytical results. In studying the potential resulting from a 
MSSM D-flat direction with supergravity induced SUSY breaking we note that at least 
for the studied parameter values, the stabihty of Q-balls requires that Q ~ 10^"^. This 
quite a large number shows that one should be careful to check whether the Q-balls 
formed in the early universe will survive scalar decays for a given set of parameters if 
Q-balls are to be cosmologically significant. 

The evaporation of Q-balls was considered for a thin-walled Q-ball and for realistic 
Q-ball profiles. We have recreated a previous result calculated for a semi-infinite 
ball. For realistic Q-ball profiles it was shown that as Q-balls lose their charge due 
to evaporation they evaporate at a faster rate. Evaporation then progresses at an 
accelerating rate until scalar decays begin. The emitted power also increases with 
decreasing charge so that Q-balls disappear in a burst of particles. Cosmologically this 
means that if the universe was dominated by Q-balls at some point in its evolution, 
there will be some inhomogeneity introduced by the evaporation process. For example, 
if the baryon asymmetry or dark matter was stored in Q-balls, inhomogeneity will be 
introduced as the Q-balls evaporate. This inhomogeneity could be cosmologically 
significant and possibly observable. 
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